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Abstract 

In this paper we show that the supergravity equations describing both cosmic 



> 

in 
10 

^ bilhards and a large class of black-holes are, generically, both Liouville integrable 

^ as a consequence of the same universal mechanism. This latter is provided by the 

O Liouville integrable Poissonian structure existing on the dual Borel algebra Bj^ of 

the simple Lie algebra An-i- As a by product we derive the explicit integration 
• . algorithm associated with all symmetric spaces U/H* relevant to the description of 

. ^ time-like and space-like p-branes. The most important consequence of our approach 

^ is the explicit construction of a complete set of conserved involutive hamiltonians 

{t)a} that are responsible for integrability and provide a new tool to classify flows 
and orbits. Wc believe that these will prove a very important new tool in the 
analysis of supergravity black holes and billiards. 



t This work is supported in part by the Italian Ministry of University (MIUR) under contracts PRIN 

2007-024045. Furthermore the work of A.S. was partially supported by the RFBR Grants No. 09-02- 
12417-ofi_m , 09-02-00725-a, 09-02-91349-NNIO_a; DFG grant No 436 RUS/113/669, and the Heisenberg- 
Landau Program. 



Contents 



1 Introduction [2] 

2 A new view-point from old results |5] 

2.1 The Ptolemaic system: Nomizu connection on a normal metric solvable Lie 
algebra H] 

2.2 The Copernican Revolution: Poissonian structure of S and Liouville inte- 
grability [8] 

2.2.1 Liouville integrability M 

2.2.2 Scenario M 

2.2.3 Realization of this scenario on B(ylN-i) [ID 

2.2.4 Involutive hamiltonians [13] 

2.3 Integrability and the metric < , >p^g [14] 

2.4 Integrability and solvable subalgebras 

3 Triangular embedding in SL(N; ]R)/SO(p, N — p; M) and integrability of the 
Lorentzian cosets U/H* \TE\ 
3.1 The integration algorithm for the Lax Equation [20] 

3.1.1 Spectral types [20] 

3.1.2 The Kodama integration algorithm for so^p n^^^"'-r) i^svisited [22] 

4 The paradigmatic example: SL(3; M)/S0(1, 2; M) versus SL(3; M)/S0(3; M) [21 

4.1 so(l, 2; M) decomposition of the SL(3; M) Lie algebra [25] 

4.2 The Lorentzian Lax operator [23 

4.3 The hamiltonian functions in involution [2H] 

4.4 Normal form of the Lax operator [29] 

4.4.1 Spectral type k = M 

4.4.2 Spectral type k = l [30] 

4.5 Characterization of orbits through the hamiltonians [31] 

4.5.1 The hamiltonians in spectral type k = [32] 

4.5.2 The hamiltonians in the spectral type k = 1 [33] 

5 Examples of explicit solutions 155] 

5.1 An example of solution of spectral type k = [33] 

5.2 Another solution of the spectral type k = with finite hamiltonians .... [36] 

5.3 An example of spectral type k = 1 [39] 

6 Conclusions [41] 



1 



1 Introduction 



Explicit supergravity solutions of pure and matter-coupled supergravity in diverse dimen- 
sions play an important role in the study of solitonic and instantonic states of superstring 
theory, in particular p-brane states [U |2] . 

Indeed one large, diversified and important class of supergravity solution^ is provided 
by the p-brane ones that are divided in two subclasses: 

• The space-like p-brane solutions that have an Euclidian world- volume and are time- 
dependent, all fields being functions of the time parameter t. 

• The time-like p-brane solutions that have a Minkowskian world volume and are 
stationary, the fields depending on another parameter t, typically measuring the 
distance from the brane. 

A view-point independently introduced in and [5j , and systematically developed in [6] , 
[Zl, [H], [9], [ID], identified, at least for space-branes, the field equations of supergravity 
corresponding to such solutions with the geodesic equations on the corresponding moduli 
space that is mostly a homogeneous space and most frequently also a symmetric space 
U/H. This identification allowed the in-depth study of supergravity cosmic billiards [TTj . 
[12j . [13], [H], [15] and lead to the discovery of their complete integrability [9], [TO] . 

It was an idea already circulating for some time in the community that also the 
construction of time-like p-branes, in particular rotational symmetric black-hole solutions, 
could be reduced to the problem of geodesic motion on appropriate moduli spaces that 
would, this time, be Lorentzian rather Euclidean coset manifolds U/H*. This idea found 
a precise formulation in the recent publication [16] . 

In connection with these applications, the question of integrability of the differential 
systems of equations describing geodesic motion on homogeneous spaces and in particular 
on symmetric non-compact cosets U/H, acquires particular relevance. As we emphasize 
in section |2| this question is intimately related with the issue of normed solvable Lie 
algebras, namely solvable Lie algebras S equipped with a non-degenerate norm < , >, 
which is positive definite in the space-brane (=billiard) case and indefinite in the time- 
brane (=black hole) case. 

In this paper, by performing a change of logical reference frame that replaces the 
route from geometry to Lie algebra into the opposite one and by gluing together pieces 
of mathematical knowledge dispersed in the literature, we show that: 

1. The integrability of all the various homogeneous models, both Euclidian and Lorentzian 
follows from the Liouville integrability of a universal parent model, associated with 
the Borel subalgebra Bn of the An-i Lie algebra. The integrability of the parent 
extends to its children algebras S if the always existing embedding S ^ Bn is 
adequate. 

^For a review and for a large set of references see for instance 
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2. Liouville integrability of Bn is an intrinsic property of this algebra which allows to 
construct an adequate number of universal hamiltonians {i)a} in involution. 

3. The norm < , > on any solvable Lie algebra S is not an independent external datum, 
rather it is intrinsically defined by the restriction to S of the unique quadratic 
hamiltonian f}o on Bn, once the embedding 5 Bn has been defined. 

4. All symmetric coset models U/H* defined as follows have integrable geodesic equa- 
tions. The Lie algebra U of the numerator is non-compact and the Lie algebra M* 
of the denominator is any of the real sections contained in U of the complexification 
He of EI C U, the former being the maximal compact subalgebra of the latter. 

5. The explicit integration algorithm has a universal form. 

Based on our new view-point we also present a new algorithmic approach to the study of 
the (eventual) integrability of homogeneous normal spaces that are not symmetric spaces, 
leaving however the actual use of such an algorithm to future publications. 

2 A new view-point from old results 

In this section we first summarize the basic facts about the Riemannian or pseudo- 
Riemannian structures that can be defined on a normed solvable Lie algebra. Our goal is 
that of reviewing the construction of the so named Nomizu connection and of its associated 
geodesic differential equations. The reason is that we aim at a Copernican Revolution. 
In this context, the classical route was from Riemannian geometry to Lie algebra theory 
since the problems that motivated the consideration of such mathematical structures were 
differential geometric in nature: in particular the geometry of scalar manifolds appearing 
in supergravity theories. It was very helpful and rewarding to find a translation vocabu- 
lary that allowed the reformulation of Riemannian geometry into a purely Lie algebraic 
setup. Yet, in relation with integrability, this classical route obscures one relevant fact: 
integrability (when it exists) is an a priori intrinsic property of the solvable Lie algebra. 
This property is intelligently, yet secretly, utilized by the (pseudo)-Riemannian structures. 
Hence, following our announced Copernican Revolution, we aim at reverting the route, 
going from solvable Lie algebra theory to (pseudo)-Riemanian geometry, rather than vice- 
versa. This change of reference frame will prove very helpful in view of old mathematical 
results, that were a little bit known in the literature on non-linear science [17], |T8], |19] 
but which had so far completely escaped consideration in the current supergravity and 
superstring literature. 

In the next subsection we prepare our Copernican Revolution with a short review of 
the Ptolemaic system. 
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2.1 The Ptolemaic system: Nomizu connection on a normal 
metric solvable Lie algebra 

Let us consider a solvable Lie algebra S. For instance S can be the Borel subalgebra of 
a complex semi- simple Lie algebra Gc, nameljj^ 

S = B(Gc) = span {Hi, (2.1) 

or it can be the solvable Lie algebra canonically associated with the pair made by a real 
form Gk of Gq and by its maximal compact subalgebra He C Gk 

S = Solv(GM/e,) . (2.2) 

Other relevant choices of the solvable Lie algebra S can be made among those associated 
with the classification of homogeneous special geometries that appear in the coupling to 
matter of supergravity theories with eight supercharges in D = 5, D = 4 and D = 3 
dimension^ 

S such that (exp[iS] , < , >) = SUGRA special Riemannian manifolds . (2.3) 

The above writing refers to the main point of the Ptolemaic system namely to the notion 
of normed metric solvable Lie Algebras. Following the original viewpoint of Alekseevsky 
we say that a Riemannian manifold {M.,g) is normal if it admits a completely solvable 
Lie group exp[Solv^] of isometrics that acts on the manifold in a simply transitive man- 
ner (i.e. for every 2 points in the manifold there is one and only one group element 
connecting them). The group exp[Solvx] is then generated by a so-called normal metric 
Lie algebra, that is a completely solvable Lie algebra Solv^ endowed with an Euclidean, 
positive definite, symmetric form < , > . The main tool to classify and study the normal 
homogeneous spaces is provided by the theorem [33], [51] that states that if a Riemannian 

^We recall that given the Cartan-Weyl basis of a complex simple Lie algebra Gc, its Borel subalgebra 
B(Gc) is defined as the solvable algebra spanned by all the Cartan generators Hi and by all the step 
operators i?" associated with all positive roots a > 0. 

■^We recall that the systematic construction of the solvable Lie algebras associated with non-compact 
symmetric spaces, pioneered in 120] and then extensively developed in the literature, has played a very 
important role in addressing, solving and systematizing a large number of supergravity problems asso- 
ciated with black-hole solutions [21], [22], [23], [24], with supergravity gaugings [25], [26] and later also 
with the issue of cosmic billiards introduced in [TT] , [T^] , [13] , [S] , [IS] and developed with the systematic 
help of the solvable Lie algebra representation of supergravity scalar manifolds in [5], [H], [S], [7], [S], [TO] . 

'^We recall that the classification of special homogeneous manifolds began with the mathematical work 
of Alekseveesky in 1975 who posed himself the problem of constructing all quaternionic Kahler manifolds 
with a transitive solvable group of isometries [2S] and then was completed and inserted into the c-map 
framework [29 of supergravity with the work of de Wit et. al. in [3D]. Further studies continued in [3T] 
and for a complete recent discussion of the topic and for all relevant further references we refer the reader 

to [3g. 
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manifold {J^,g) admits a transitive normal solvable group of isometries exp[Solv;vt], then 
it is metrically equivalent to this solvable group manifold 



M ~ exp [SolVyy4 ] , 

g\eeM = <,> (2.4) 

where <, > is the Euclidean metric defined on the normal solvable Lie algebra Solv^j. 

The conjecture of Alekseevsky was just restricted to quaternionic Kahler manifolds 
and implied that any such manifold Ai that was also homogeneous and of negative Ricci 
curvature should be normal, in the sense over mentioned, namely a transitive solvable 
group of isometries exp [So1v;k] should exist, that could be identified with the manifold 
itself. Note that the actual group of isometries U of could be much larger than the 
solvable group, 

U D exp [Solv^^] , (2.5) 
as it is for instance the case for all symmetric spaces 

M = ^ (2.6) 

yet the solvable normed Lie algebra (So1v_a4 , < , >) had to exist. The problem of clas- 
sifying the considered manifolds was turned in this way into the problem of classifying 
the normal metric solvable Lie algebras {S, < , >). Note that in Alekseevsky's case the 
symmetric form < , > was not only required to be positive definite but also quaternionic 
Kahler. Alekseevsky's conjecture actually applies to more general homogeneous Rieman- 
nian manifolds than the quaternionic ones: for instance it applies to all those endowed 
with a special Kahler geometry or with a real special one as the classification of de Wit 
et. al. [30] demonstrated. It also applies to the symmetric spaces appearing in the scalar 
sector of extended supergravities with more than eight supercharges. For all these mani- 
folds there exists the corresponding normal metric algebra {S, < , >), in other words they 
are normal. This happens because they are Einstein manifolds of negative Ricci curvature 
and, although we are not aware of any formal mathematical statement in this direction, 
one might make the 

Conjecture 2.1 << Every homogeneous Einstein manifold M. of negative Ricci curva- 
ture is normal, namely there exists a normal metric solvable Lie algebra (iS, < , >) such 



that identifying S with So\v_\4 eq. (2.4) applies. » 



Proving such a conjecture amounts to proving that for every homogeneous Einstein 
manifold the group of isometries U, which by hypothesis of homogeneity exists and 
has a transitive action on the manifold, admits a solvable simply transitive subgroujJ^ 
exp [5] C U. If this is true, in view of the already mentioned theorem the rest follows. 



""Simply transitive means that each group element has no fixed points. 
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The key assumption is the negative Ricci curvature. Manifolds of positive Ricci curva- 
ture, which are typically compact, are excluded. All compact symmetric spaces are indeed 
counterexamples. For U/H compact there is no transitive solvable subgroup of U. 

The recollection of these well known facts was done in order to emphasize the following 
point. In the Ptolemaic system that starts from Riemannian geometry and arrives at 
solvable Lie algebras S, this latter emerges in conjunction with a well defined metric form 
< , > defined over it. For instance if we focus on Borel solvable algebras B(Gc), they are 
endowed with the following canonical metric: 

< Hi , Hj > = 2 5ij , 

<ni,E''> = , 

<E^E^> = 6^f, (2.7) 

whose normalization is absolute if the generators {Hi, -E", E~°'} of the Weyl-Cartan basis 
for Gc have the standard normalization 

[Hi , Hj] = , 

[E^'^E^] = N{a,(3)E''+^ if a + /? is a root ; \N{a,(3)\ = l, 
[E" , E""] = 2a ■ n, 

rank 

a ■ (3 = caPi ^ {Hi} = orthonormalized basis . (2.8) 

i=l 



The metric (2.7) is singled out by the relation of the Borel algebra B(Gc) with one specific 



Riemannian Einstein manifold of negative Ricci curvature. This latter is 

M = ^ (2.9) 

where G^s denotes the group generated by the unique real section of the complex 
Lie algebra Gc which is maximally split (or, equivalently maximally non-compact) and 
He denotes the unique maximally compact subgroup of Gms- It turns out that B(Gc) is 
just the solvable Lie algebra of this Riemannian manifold Ai 

B(Gc) = SoIvm (2.10) 



and the metric < , > (2.7) on the solvable Lie algebra Solv^ is that induced by the unique 



Einstein Riemannian metric on the corresponding coset *^ 



He 



2.9). 



Similarly it happens in all other constructions of the Ptolemaic system. The re- 
quirements imposed on the final Riemannian Einstein manifold one wants to construct 
predetermine the metric < , > on the solvable Lie algebra. 

Once the metric form is given, the construction of geometry and of the associated 
geodesic equations follow uniquely. The issue is just that of calculating the Levi-Civita 
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connection of the metric g induced on the manifold by the form < , > defined on the 
solvable Lie algebra. One way of describing this Levi-Civita connection is by means of 
the so called Nomizu operator acting on S. The latter is defined as follows: 

L : 5 5^5, (2.11) 
VX,y,ZG5 : 2 <l.xY,Z >=<[X,Y],Z > - < X,[Y,Z]> - <Y,[X,Z]> . 

The Riemann curvature operator on S can be expressed as 

Riem(X, Y) = [Lx, W] - h[x,Y] ■ (2-12) 

If we introduce a basis of generators {T^} for S and the corresponding structure constants 
defined by 

[Ta,Tb] = Jab'^Tc (2.13) 

together with the metric tensor: 

<TA,TB>=gAB (2.14) 



the connection defined by eq.(2.11) leads to the following connection coefficients: 

^aTb = ^AB 

r^B = (2.15) 

which are constant numbers. 

Equivalently, we can define the Levi-Civita-Nomizu connection starting from the dual 
description of the solvable Lie algebra in terms of Maurer-Cartan equations. Let be a 
basis of Maurer-Cartan forms dual to the generators Tg, namely e^iTs) = 6^. We have 

de"^ = IfAB^'e^ Ae^ (2.16) 

Interpreting as the vielbein over the solvable group manifold we write the vanishing 
torsion equation 

= cie^ + A qbc (2.17) 

where = —u^^ is the standard so(?2)-Lie algebra valued spin connection [n = 

dim(iS) ). The relation between the two descriptions is immediate: 

c.^^ = r^^^^^e^ (2.18) 

where the tensor g^^ is automatically antisymmetric in force of its definition. 

Given the connection coefficients the differential geodesic equations can be immedi- 
ately written. In the chosen basis the tangent vector to the geodesic is described by 
n fields Y^{t) which depend on the affine parameter t along the curve. The geodesic 
equation is given by the following first order differential system: 

+ r^^y^F^ = . (2.19) 
The above equation contains two data: 
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1) the structure constants of the solvable Lie algebra /^^^ 

2) the metric tensor qab- 

As we emphasized the second datum, namely the metric, comes down, in the Ptolemaic 



system from the geometric interpretation of the differential system (2.19) as geodesic 
equations on a Riemannian Einstein manifold and it is the metric of that manifold what 
eventually predetermines qab- 

Let us now implement our Copernican Revolution and let us forget for a moment 



about the Riemannian structure. The system (2.19) is just a non linear differential system 
defined over the dual S* of a Solvable Lie algebra S. How could we directly derive such 
a differential system and may be established its Liouville integrability from the very 



structure of SI Is there a way of deriving the metric tensor implicitly contained in (2.19) 
from the Lie algebra SI There is. 

2.2 The Copernican Revolution: Poissonian structure of S and 
Liouville integrability 

Given the solvable Lie algebra S we can consider the co-adjoint orbits of the corresponding 
solvable group exp[iS]. In simple language this amounts to consider functions defined over 
the dual Lie algebra S* . An element of S* is a linear functional 

C : 5 ^ C or M (2.20) 

and using a general theorem in linear algebra the dual of a finite dimensional vector space 
is isomorphic to the same space. In practice, given a basis Ta of S we immediately obtain 
the dual basis by defining 

C^{Tn) = 5^ . (2.21) 

Any dual Lie algebra element 2) G 5* can be written as a linear combination of the dual 
generators 

^ = YaC^ (2.22) 

and the most general function (j) on S* is actually a function of the n coordinates Ya- 
Given any two such functions (pi and 02 we can define their Lie-Poisson bracket in the 
following manner: 

{01, 02} = fAB yCg^^g^. (2.23) 

In this way the space of co-adjoint orbits becomes a Poisson manifold, independently 
from the existence of any metric < , > on S. Then one is allowed to consider evolution 
equations of the following form: 

^Ya + {Ya, H} = (2.24) 
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where Ti = HiY) is some function on the dual Lie algebra S* that we can regard as the 
Hamiltonian. 



The question is whether the geodesic equations (2.19) can be put in the hamiltonian 



form (2.24), namely, whether there exists a hamiltonian function, necessarily quadratic, 
which reproduces the Levi-Civita-Nomizu connection. The answer is obviously yes. It 
suffices to write 

n = ^YAYBg'''' (2.25) 
and identify the variables and Yb through the relation 

Ya = 9abY''. (2.26) 



Both in eq. (2.25) and in eq.(2.26) there appears the metric qab defined by the non 
degenerate normal form < , >. Hence it may seem that we made no real progress and we 
simply rewrote the same equations in a different style. The notion of the metric tensor is 
still essential and it looks external to the pure Lie algebraic structure. It is not so, as it 
appears from the following argument. 

2.2.1 Liouville integrability 

The key point we would like to emphasize is that the definition of the Lie-Poisson bracket 



(2.23) depends only on the structure constants of the algebra S and nothing else, so it is 



intrinsic to the algebra. Let us next recall the notion of Liouville integrability. 

Definition 2.1 A symplectic manifold of dimension 2m endowed with a Lie-Poisson 
bracket {, } is Liouville integrable if there exists m functionally independent functions 
^iiy) of its 2m coordinates 3^/ that are in involution. Namely we must have: 

V2,j = l,...m : {$,,$^.} = (2.27) 

and 

rank I 

dyj 

When these conditions are fulfilled any of the functions $j can be chosen as the hamilto- 
nian Ti. and the corresponding Euler equations 

= + {yi,n} (2.29) 

at 

are completely integrable, since they admit m first integrals of the motion $i(3^). 
Let us now envisage the following: 



rank — f ) = m . (2.28) 
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2.2.2 Scenario 

Given the solvable Lie algebra S, whose dimension we denote 

ds = diniiS , 



(2.30) 



imagine that with respect to its intrinsic Lie-Poissonian structure (2.23) there exist ps 



functions i)aiY) of the coordinates Ya {A = 1, ...,ds) on S* that are in involution 

{f),, i)fs} = (a,/3 = l,...,ps) (2.31) 
with the integer ps lying in the range: 



ds > Ps > 



d^ 
2 



and suppose furthermore that, having defined the integer 

c = 2ps - ds , 



(2.32) 



(2.33) 



precisely c of the functions i)a{Y), which, for this reason we rename Ce(Y) (£ = 1, . . . c), 
are Casimirs, namely they commute with all the coordinates: 

{YA,Ci} = ; (i = l,...,c; A = I,... ds) . (2.34) 

Under these hypotheses, by rearranging the set of involutive functions as it follows: 



^i=l,...,ms e=l,...,c. 



(2.35) 



we obtain the following situation. The level surfaces ^ri,...,rc defined by setting the c 
Casimirs to fixed values 



CM 



[ CM 

are by definition manifolds of even dimension: 

dim^ri,...,rc = 2m5 

ds - c 
ms = — - — 



ri E R (or C) 
r2 G M (or C) 

G M (or C) 



ds - Ps 



(2.36) 



(2.37) 



and on these manifolds there exist exactly ms functions in involution, namely the ^i{Y) 



of eq.(2.35). It follows that each of these manifolds acquires a symplectic hamiltonian 
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structure. Naming yj the 2 ms coordinates on ^ri,...,rc and choosing as hamiltonian any 
hnear combination of the pull-backs on ^ri,...,ra of the functions ^i{Y) 

n{y) = pull-back on %.,,...,r, of a' (2.38) 

the corresponding Euler equation differential system 

^ j^yi + {yi,n} (2.39) 

is by construction Liouvillc intcgrable. Indeed the pull-backs of the ms functions ^i{Y) 
provide the necessary first integrals of motion. 

2.2.3 Realization of this scenario on ]B(An-^i) 

The above described scenario is precisely realized in a very notable case, namely that of 
the Borel subalgebra B(ylN~i) C A-^-i- As it is well known the simple Lie algebra An-i, 
identified by the following Dynkin diagram: 

o — o — 0---0 — o — o 

«! a2 tta OiN-S OiN-2 OiN-l 



is the abstract form of the Lie algebra s[(N; C) of complex traceless matrices in dimension 
N. The corresponding Borel subalgebra Bn = B(74n_i) is simply given by the subset of 
all upper triangular traceless matrices. It is more convenient to relax the condition on 
the trace and consider the Borel subalgebra Bn = B (0((N; C)) which is simply made by 
all upper triangular matrices. Reduction to Bn will be performed by putting one of the 
Casimirs to the null value. Hence we define 



0[(N;C) D 



3 b 



( " 


-k 


-k 


-k 


■ " 1 





-k 


-k 


-k 


-k 








-k 


-k 


-k 











-k 


-k 



y 

The dimension of this solvable algebra is easily computed 

N(N + 1) 



dim] 



"N 



(2.40) 



(2.41) 
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N 


c/n 


Pn 


cn 


2m 


2 


3 


2 


1 


2 


3 


6 


4 


2 


4 


4 


10 


6 


2 


8 


5 


15 


9 


3 


12 


6 


21 


12 


3 


18 


7 


28 


16 


4 


24 


8 


36 


20 


4 


32 


9 


45 


25 


5 


40 


10 


55 


30 


5 


50 


11 


66 


36 


6 


60 


12 


78 


42 


6 


72 


13 


91 


49 


7 


84 


14 


105 


56 


7 


98 


15 


120 


64 


8 


112 


16 


136 


72 


8 


128 



Table 1: In this table we give the dimensions of the Borel algebras B(g[(N)), the corre- 
sponding number of functions in involutions, the corresponding number of Casimirs and 
the ensuing dimensions of the orbits for the first sixteen values o/N. 

It was demonstrated by Arhangel'skii in [17J that this Poissonian manifold is Liouville 
integrable according to the scheme described in the previous subsection and the procedure 
to construct the required set of functions in involution was described in [TTl [18] . Let us 
review these results. We distinguish the two cases where N = 2i/ is even and where 
N = 2i/ + 1 is odd. The number of functions in involution, the number of Casimirs 
and the ensuing even dimension of the orbits is displayed below: 





d = dimi 


p = T^funct. in inv. 


c = T^of Casim. 


2m = dim of orbits 




v{2v + l) 


u"^ + 1/ 


V 




^2u+l 


{u+l){2v + l) 


v"^ +2u+l 


y + l 


2{iy^ + u) 



(2.42) 

For the reader's convenience in Tabled] we tabulated the first instances of such numbers. 
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2.2.4 Involutive hamiltonians 



How are the functions i)a(X) explicitly constructed? 

First let us observe that a simple set of coordinates on the Borel subalgebra B(0[(N)) 
is simply given by the entries of the upper triangular matrix b mentioned in eq.( |2.40[ ). 
Secondly, for reasons that will become clear later on, let us consider an N x N matrix L 
which satisfies the condition 

(L 7]f = Lr] (2.43) 

where 

r/ = diag(-l,+l,... ,-!,+!, +1,... + 1) , p<q ; p + q = N (2.44) 
is any of the available choices of {p, q) signatures in dimension N. Irrespectively of the 



choice of ?], the number of parameters contained in a matrix L satisfying (2.43) is always 
equal to N(N + l)/2 which is the dimension of the Borel algebra B(g[(N)). Indeed there 
is a simple one-to-one map from the space of upper triangular matrices to the space of 



matrices satisfying eq.(2.43), which reads as follows 



b 1-^ L = b + rfh i] , 

L ^ b = L> - ^ diagL . (2.45) 

In the above equation we have used the following convention. For any matrix M, we 
denote by M> its upper triangular part and by M< its lower triangular one including the 
diagonals. 



Hence the space of matrices fulfilling eq.(2.43) provides a coordinate basis for the Borel 



algebra. Thanks to the simple relation (2.45), the use of different i] tensors just amounts 



to a linear coordinate transformation on Bn, in other words to a change of basis for the 
generators of the solvable Lie algebra. What is then the relevance of considering such 
different choices of r] rather than focusing on a single conventional one? The answer is 
elementary. Each different rj prepares a basis well adapted to the decomposition of the Lie 
algebra 0t(N) with respect to its subalgebra so{p, q). Such a decomposition is necessary in 
order to study the geometry and the associated geodesic equations of the coset manifold 

GL(p + q) 

All these (p, q) systems are integrable since all the manifolds M.p,q are metrically equivalent 
to the same solvable manifold exp[Bp+q], equipped however with a different normal form 
< , >p,g, which is positive definite only for {p = 0, q = N}. The Liouville integrability of 
all these systems has the same common root, namely the existence of the p^ independent 
functions in involution that now we display and which is an intrinsic property of the Borel 
algebra Bp+g. 
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The algorithm of constructing these functions, originally derived for {p = 0, g = N} 
in [17] and [18], being generalized to the case under consideration, is the f ollow ing. 

Starting from the parameterization of Bn by means of the matrix L (2.45) fulfilling 



eq.(2.43) the complete set of pn functions [)a that are involutive with respect to the 



Lie-Poisson bracket (2.23) is enumerated by an ordered pair of indices 

(a, 6) 



a 



where: 



a 
b 







'N" 
_2'_ 
1 N-2a 



The functions f)^^ can be iteratively derived from the following relation: 

det {{L - X)ij : a + l<z<N, l<j<N-a} 

N-2a 

^ao I + ^ t)., A^-^'^-" 1 , a = 0,... 



b=l 



(2.47) 



(2.48) 



(2.49) 



where, by definition, Sao is the coefficient of the power A^"^*^. 

As we know from the previous discussions cn of these generalized hamiltonians are 
actually Casimirs. 

For example, for the case N = 3 that we study in detail in the next sections the 
functions \)q2, f)o3 are pure hamiltonians while f)oi, f)ii are Casimir functions. 



2.3 Integrability and the metric < , > 



Having established in purely algebraic intrinsic terms the integrability of the differential 
system based on the solvable Lie algebra S = Bn, we can now answer the question about 
the origin of the metric form < , >p ^ which, in the Ptolemaic system, leads to the very 
same differential equations. The problem is very simply solved by recalling the quadratic 



form (2.25) of the hamiltonian which yields the geodesic equations associated with the 
Levi-Civita-Nomizu connection. The metric qab is defined by the coefficients appearing 
in general unique quadratic hamiltonian function 7i = f)o2- This view-point is actually 
very efficient. Changing the metric signature rj amounts, as we have emphasized, to a 



change of basis on the solvable Lie algebra. Using all the time the same formula (2.49) 
to calculate the hamiltonians ()a, but varying the choice of t] in the definition of L (2.45), 
results in a change of the coefficients in the quadratic hamiltonian \jao 5 namely in a change 
of metric qab over the Lie algebra. All the various coset models with indefinite signature 
are covered in this way by a unique algorithm. 
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2.4 Integrability and solvable subalgebras 

In view of the above results on the integrabihty of all Borel algebras Bn there arises a 
natural question about the relation of this integrability with the possible integrability of 
its subalgebras. Indeed in view of the theorem that states that every linear representation 
of a solvable algebra S admits a basis where all the elements X & S are represented by 
an upper triangular matrix it follows that any solvable Lie algebra can be regarded as a 
subalgebra S C Bn for a suitable choice of N. ^ 

Hence let us consider a generic subalgebra S C Bn and name t)a the pull-back on S* 
of the pn generalized hamiltonians defined on Bn- In particular we are just interested in 
the differential system constructed on S* where the Hamiltonian Ti is the pull-back of the 
unique quadratic hamiltonian Ti = f}o2 constructed on Bn- From the mathematical point 
of view there might be other choices but from the physical point of view this is the only 
relevant one. So let us put 

n = 1)02 (2.50) 
and let us divide the Bn coordinates in two subsets, those along S and those normal to S 

{Ya} = W^} . (2.51) 

Next let us compare the evolution equations of the subalgebra S with the pull-back on S 
of the Btv equations. In view of our choice (2.50) we have 

n{x) = nix,w)\w=o = n{x,o) . (2.52) 

The 5-equations are simply read 

jx, = -{x,,n} ^ 

= -f.^XkOm. (2.53) 
while the pull-back on S of the Bn equations is 

~T:^i = —{Xi,'H}\w=o 
at 

= -f.^'Xkdm - fJXkd''n\w=o , (2.54) 

= -faf3'Xkd''n\w=o - fJXkd'n . (2.55) 

The conditions necessary for the system ( |2.54|2.55 ) to reduce consistently to the system 
(2.53) are the following ones: 

d^n\w=o = 0, (2.56) 
fa' = 0. (2.57) 
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The first condition is satisfied if there are no mixed coefficients in the metric imphcitly 
defined by the hamiltonian, namely if 



9' 







(2.58) 



The second condition implies that the decomposition of the Borel algebra Bn with re- 
spect to its solvable subalgebra S should be reductive. Altogether the conditions for 
the consistent reduction of the system ( 2.54|2.55 ) to the system (2.53) can be written as 
follows: 



Bn = 5 ( 
<S ,S^> = , 
[S,S] c 5, 



5, 



(2.59) 
(2.60) 
(2.61) 
(2.62) 



Under the hypotheses (2.59-2.62) a geodesic on the manifold (exp[S'], <, >) satisfies 
the geodesic equations on the embedding manifold (exp^^j, < , >). The latter is a 
symmetric space, namely GL(N, M)/SO(p, N — p), so that (exp[S'], < , >) is revealed to 
be a geodesically complete submanifold of a symmetric space. In force of a general theorem 
[M] this implies that also (exp[S'], < , >) is a symmetric space. 

Hence the exceptional cases of homogeneous solvable manifolds (exp[5], < , >) that 
are not symmetric spaces (compare with the classification of homogeneous special geome- 
tries, |30]) are based on solvable Lie algebras S whose embedding in Bn is not reductive 



and violates eq.s (2.59-2.62). In order to establish their Liouville integrability one has 



to study the pull-back on S of the pn involutive hamiltonians (see Table ij) and as- 
certain how many of them remain in involution among themselves and with the unique 
hamiltonian H = i)o2- Such a study involves several complicacies and it is postponed to a 
subsequent pubhcation. For the case of symmetric spaces, instead, the proper reduction of 
the hamiltonians and the consequent integrability is guaranteed a priori from the proper 
reduction of the Lax representation that we discuss in the next section. 



3 Triangular embedding in SL(N; M) /SO(p, N - p; R) and 
integrability of the Lorentzian cosets U/H* 

As already recalled in the introduction it was realized in recent years that the field equa- 
tions of supergravity describing either 

1) space-like p-branes as the cosmic billiards, or 

2) time-like p-branes as several rotational invariant black-holes in = 4 and more 
general solitonic branes in diverse dimensions 
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reduce to geodesic equations on coset manifolds of the type 

M^l or M*^^^ (3.1) 

where U = exp[U] is the group manifold generated by the duality algebra U relevant 
to the considered supergravity model in the considered dimensions. As a rule without 
exceptions U is always some non- compact real form of a complex Lie algebra Uc- The 
Lie algebra EI of the subgroup H is instead the unique maximal compact subalgcbra of U. 
The coset A4 corresponds to the case of space-like p-branes, in particular cosmic billiards. 
In the second coset M.*, pertaining to the case of time-like branes, the Lie algebra H* of 
H* is another non-compact real section of the complexification He, while U remains the 
same. It happens indeed the following situation always occurs: the real duality algebra 
U admits as proper subalgebras a few different instances of real sections of Mq which 
necessarily include the maximally compact one H. 

We recall few illustrative examples. Consider the case of maximal supersymmetry 
with 32 supercharges. The duality algebra in D = 4 is U = £7(7) whose maximal compact 
subalgebra is H = su(8) C £7(7). The complexification of this is He = 5((8,C). Aside 
from 5u(8) the other real sections of sl(8, C) which are subalgebras of £7(7) are the following 
ones: = 5t(8,M) C £7(7) and H* = su'^(8) C £7(7). For the same number of 
supercharges the duality algebra in D = 5 is U = £5(6) • Here we have H = U5p(8), 
He = 5p(8,C), H* = usp(4,4) or H* = sp(8,R). 

Consider instead the case of 1/2 maximal supersymmetry, namely supergravity with 
16 supercharges. In D = 3 the theory with 6 + n supervector multiplets has the following 
duality algebra U = so(8,8 + n). The corresponding maximal compact subalgebra is 
H = so (8) X so(8 + n) and the available H* are: H* = so{8 — p,p) x so{8 + n — p,p) 
(p=l,2,3,4). 

For all these cases we can make the following statement. 

Statement 3.1 << Let N be the real dimension of the fundamental representation ofU. 
For each choice o/H or H* there exist a suitable integer p < [^] and a diagonal metric 

T) = diag( ^-l,+l,. . .,-l,+l , +1,+1^,... + 1 ) , (3.2) 

2p N-2p 

such that we have a canonical embedding 

U ^ s[(N;M) , 
U D H* 50(p,N-p;M) C s[(N;M) . (3.3) 

This embedding is determined by the choice of the basis where Solv (U/H*) is made by 
upper triangular matrices. In the same basis the elements ofK are rj-symmetric matrices 
while those o/H* are r]- antisymmetric ones, namely: 

VK G K : r]K^ = K^r] , 

VH e H* : r/H^ = -H^ry . (3.4) 
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Just as in the case of cosmic billiards the embedding of the Riemannian coset U/H into 
the universal covering coset SL(N; ]R)/SO(N; M) provided the key to obtain an explicit 
integration algorithm for the associated first order geodesic equations, in the same way the 



embedding (3.3 ) of the pseudo- Riemannian U/H* into SL(N; M) /SO(p, N — p; M) provides 
the key to extend the same integration algorithm also to this indefinite metric symmetric 
spaces. Indeed that algorithm is defined for SL(N; ]R)/SO(p, N — p; M) and it has the 
property that if initial data are defined in a submanifold U/H* where U C SL(N;]R) 
and H* C SO(p, N — p; M), then the entire time flow occurs in the same submanifold. 
This property is a simple consequence of the following argument. Let L G K be an 77- 



symmetric matrix satisfying property (2.43) which lies in the subspace K corresponding 



to the symmetric decomposition U = H* © K of the Lie subalgebra U C g[(N;]R). By 
construction the corresponding connection W = L> — L< lies in H*, namely we have 
W e B*. It follows that [W, L] e K. Hence if L(0) G K, the evolution equation 
L + [W, L] = never brings L(t) out of that space. 

Hence the embedding ( 3.3[ ) suffices to define explicit integration formulae for all su- 



pergravity time-like p-branes based on pseudo-Riemannian symmetric spaces. 
Let us review the steps of the procedure. 

1. First one defines a coset representative for U/H* just as in the Riemannian case, 
namely: 

/=i 

L(0) = Yl ew[<^iE"'] exp[hiW] (3.5) 

I=m 

where the roots pertaining to the solvable Lie algebra are ordered in ascending order 
of height (a/ < aj if I < J), W denote the non compact Cartan generators and 



the product of matrix exponentials appearing in (3.5) goes from the highest on the 
left, to lowest root on the right. In this way the parameters {0} = {(fj,hi} have 
a precise and uniquely defined correspondence with the fields of supergravity by 



means of dimensional oxidation [5l |BJ. Thanks to the mapping (2.45), the same 
upper triangular matrix L ((f)) can be regarded as a coset representative for U/H or 
for U/H*. 

2. Restricting all the fields (p of supergravity to pure time dependency (p = 0(t), the 
coset representative becomes also a function of time L (0(t)) = h{t) and we define 
the Lax operator L{t) and the connection W{t) as follows: 



L{t) = 5^TrfL-i^LK, )K, 



W{t) = ^ Tr (^L-i^LH,^ H, (3.6) 



^In the case of time like p-branes such as black-holes time is actually the radial coordinate. 



18 



where Kj and denote an orthonormal basis of generators for K and H*, respec- 
tively 



Tr (K, , K, 



Tr(Hi,H, 
Tr(Ki, H 




n = 

m = dim (H*) = dim (H) , 

r = # compact generators of K = 7^ non-compact generators of EI . 

(3.7) 

3. With these definitions the field equations of the supergravity time-like p-brane or 
black-hole, which are just the geodesic equations for the manifold U/H* in the 
solvable parametrization, reduce to the single matrix valued Lax equation ^ 

+ [Wit) , Lit)] = (3.8) 

which is the compatibility condition for the linear system exhibiting the iso-spectral 
property of the Lax operator L 

L ^ = ^ A , (3.9) 

^^ = 1^^ (3.10) 

dt ^ ' 

where 

A = diag(Ai,...,A^) (3.11) 
is the diagonal matrix of eigenvalues and \l/(t) is the eigenmatrix. 

4. If we are able to write the general integral of the Lax equation, depending on 
n = dim(U/H*) integration constants, then comparison of the definition of the Lax 
operator ( 3.6[3.5 ) with its explicit form in the integration reduces the differential 



equations of supergravity to quadratures 
d 



(j)it) = Fit) = known function of time. (3.12) 
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3.1 The integration algorithm for the Lax Equation 



Let us assume that we have exphcitly constructed the embedding (3.3). In this case, in 
the decomposition 

U = K©e* (3.13) 

of the relevant Lie algebra U, the matrices representing the elements of K are all t]- 
symmetric while those representing the elements of M* are all r^-antisymmetric as we 
have already pointed out. Furthermore the matrices representing the solvable Lie algebra 
Solv(U/H*) are all upper triangular. These are the necessary and sufficient conditions 



to apply to the relevant Lax equation (3.8) the integration algorithm originally described 



in [19] and reviewed in [9l [TO]. The key point is that the connection W{t) appearing in 



eq.(3.8 ) is related to the Lax operator by means of the already recalled algebraic projection 



operator as follows: 

= n(L) := L> - L< . (3.14) 



The relation (3.14) is nothing else but the statement that the coset representative L(0) (= 
L + W) from which the Lax operator is extracted is taken in the solvable parametrization. 

The only relevant new feature distinguishing the indefinite metric case from the definite 
one is the discussion of the spectral types. 

3.1.1 Spectral types 



From eqs.(2.43) it follows that, generically, the Lax operator L(0) is not a symmetric 
matrix. It is such only in the Euclidean case {p = and g = N). Therefore L{0) 
eigenvalues are generically complex numbers. We will concentrate on the case when L{0) 
is a simple matrix, i.e. all its eigenvalues are distinct. Then, in order for L{0) to be real, its 
eigenvalues have to group in k complex conjugated pairs and N— 2A; real eigenvalues, where 
k is some fixed integer in the range < k < p. Obviously, Lax matrices corresponding to 
different values of k can not be related by a similarity transformation. Thus the integer k 
parameterizes inequivalent initial data Lk{0) that we decide to name spectral types. For 



the metric (2.44) one can choose a basis in the space of eigenvalues where their complex 
conjugation properties become 

^*2a-l = ^2a , a = l,...,k, 0<k <p, 

a; = A„ , a = 2fc + 1,...,N . (3.15) 



Simple inspection of eqs. (3.15) shows that the corresponding matrix eigenvalues (3.11 ) 



and its complex conjugated matrix A^, are actually related by a similarity transformation 

Afc=TfcAfcT-i (3.16) 

which will be useful in what follows. Here, the complex symmetric N x N-matrix Tk has 
a very simple block-diagonal structure 

Tfc :=diag(5,...,5,-lN-2fc) (3.17) 
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and satisfies the properties 

Tfc Tfc = In , Tkr]Tk = r]. (3.18) 
The matrix comprises k sub-blocks given by the (2 x 2)-matrix B, defined below: 



B 




B^ = -1, 



(3.19) 



and the bottom ((N — 2k) x (N — 2/c))-matrix sub-block which is proportional to the 
((N — 2k) X (N — 2/c))-unity matrix lN-2fc; i is the imaginary unity. The matrices 
{k = 0, 1, ...,p) satisfying eq. (3.16) form the initial data for the L(0) eigenvalues. 



Now, let us turn to constructing the initial data for the L(0)-eigenmatrix \E'(0) (3.9). 
If the Lax eigenvalues are complex (i.e., ^ Ak), then the eigenmatrix ^1^(0) has to be 
complex as well in order for the Lax matrix 



to be real 



Lfc(O) = ^fc(O) Afc ^^\0) 



^fc(O) Afc (^fc(O)) =^fc(0) Afc^-i(O) 



(3.20) 



(3.21) 



Besides the constraint (3.21), the complex eigenmatrix \E'(0) which diagonalizes the initial 
Lax operator L[0) should satisfy one more constraint [I9j 



^i{o)v^k{o) = v 



(3.22) 



namely it should belong to the complexified group SO(p, N — p; C). Let us introduce the 
new set of block-diagonal complex, symmetric N x N-matrices 



T, 



k ■'- 



diag 



-j= (I2 - B), -j= (I2 - B), lN-2k 



(3.23) 



that satisfy the following relations among themselves, with the metric rj (2.44) and with 
the previously introduce set (3.17): 

Tk Tk = In 



kVTk = V: TkTk = -Tk 



(3.24) 



Using eqs.(3.16) and (3.24) by straightforward calculation we can verify that the complex 
eigenmatrix 



^fc(O) := Oo Tk 



(3.25) 



satisfies both constraints (|3.21|) and (|3.22|) if the introduced real matrix Oq, entering into 

(3.26) 



eq. (3.25), satisfies the pseudo-orthogonality condition 
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namely if Oo e SO(p, N - p; M). 

We conclude that the initial data for the Lax operator are represented in the following 
way: 

Lfc(O) = Co fk Ak f^^ O^^ where 

Co e SO(p,N-p;M) . (3.27) 

Hence the manifold of solutions of the Lax equation for the case GL(N, ]R)/SO(p, N — p; M), 
splits into (p + 1) disconnected branches corresponding to the spectral types = 0, 1, 

When we consider a subsystem U/H*, having fixed p, the actual number of spectral 
types is determined by considering which normal forms 

Afc = % Afc f^^ (3.28) 

actually belong to the space K. Indeed the spectral type k will be included in the inte- 
gration algorithm if and only if 

Afc G K . (3.29) 

The above discussion of spectral types actually coincides with the discussion of normal 
forms reported in [loj^ In the same paper the normal form was group-theoretically 
interpreted as an element of the subspace IK belonging to a subvector space of the form: 

Furthermore in the number of actually available normal forms (spectral types in our 
nomenclature) was shown to admit the following group theoretical interpretation: 

# of spectral types = rank ^— ^ + 1 (3.31) 

where C H* is the maximal compact subgroup of H* . 

3.1.2 The Kodama integration algorithm for gQ^p^^-'p-R) revisited 

Having clarified the fundamental issue of spectral types let us describe in full detail the 
adaptation of the Kodama integration algorithm [19] to the indefinite metric case. 

Summarizing all our previous discussions the starting point of the algorithm is provided 
by the initial data listed below: 

a) The spectral type, codified by the choice of one of the p + 1 matrices {k = 0, 1, ...p) 
(3.23). At fixed k we set: 

r := ffc . (3.32) 



''We would like to thank our friend and long term collaborator M. Trigiante who attracted our attention 
to the constructions of paper ^16j when the results of the present paper were already in final form. 
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b) The eigenvalue vector: 



N-l 



A 



Ai, A2, . . . , Atv-i , — Aj 



(3.33) 



i=l 



If the spectral type is k — all the eigenvalues are real. If the spectral type is 
k — r then the first 2r eigenvalues are arranged in complex conjugate pairs such 
that — A2i-i while the remaining N — 2r are real. 

c) The choice of an arbitrary element of the pseudo-rotational group Oq G SO(p, N — p; M) , 
namely a real N x N matrix satisfying: 



In terms of these initial data we define the following complex matrices: 

^(0) = OqT e SO(p,N-p;C) . 



(3.34) 



(3.35) 



Next we construct a time- dependent N x N matrix c{t) whose elements are defined by the 
following formula: 



N 



Ci^ W = 5^ *ife(0)exp[-2Re(Afe)i] (cos[-2Im(Ajk) + isin[-2Im(Ajk) i]) (*-'(0))^^. , 

(3.36) 



k=l 



and we introduce the N-functions Di(t) (-Do(t) := 0) constructed as the determinants of 
the principal diagonal i x i sub- matrices of c{t), namely: 



De{t) = Det 



it)) 
/ i<i,j<e 



(3.37) 



Using these building blocks we construct a time-dependent, real, pseudo-rotation matrix 
0{t) in the following way. First define a time-dependent complex matrix ^'(i) whose 
entries are given by 



, , exp[— ReAj t] (cos[— ImAj t] + i sin[— ImAj t] ) 



and then we obtain the desired 0{t) by setting 

0{t) = ^{t)T-^ . 



^ Cll ••• ci,i_i ^'1^(0) ^ 

\cii ... Ci,i_i ^'ij(O) I 

(3.38) 



(3.39) 



23 



It is a matter of direct verification that defined as above, the evolving SO(p, N — p) group 
element 0{t) is indeed real at all instants of time 

0{t) = 0{t) , (3.40) 

thus 0{t) G SO(p,N — p;M). Finally the explicit form of the Lax operator solving Lax 
equation with the chosen set of initial conditions is given by 

L{t) = 0{t)AO-\t) , (3.41) 
A = TAT-^ (3.42) 

where: 



A 



/A, 










\ 





A2 



























Aat-I 





V 













(3.43) 



At first sight the reader might consider baroque the substitution 

exp[-2Xjt] exp[-2Re(Afc)t] (cos[-2 Im(Afe) t] + i sin[-2 Im(Afc) t]) 



(3.44) 



used both in eq. (3.36) and eq. (3.38). Actually such a substitution is very handy in order 
to verify the reality of the solution and also in order to understand its analytic structure. 
First of all it is fairly simple to check that the matrix Cij{t) and hence all of its minors are 
real. Secondly thanks to the same token one verifies that also 0{t) is real and all of its 
entries are rational functions of exponentials, cosines and sines or square-roots thereof. 
The periodic trigonometric functions are absent for the spectral type A; = and appear 
only for k > 1. The appearance of these sines and cosines is the truely new feature due 
to the pseudo-Riemannian structure of the denominator group H*. 

For all spectral types the evolution of the Lax operator is given by a time dependent 
SO(p, N — p; M) similarity transformation, starting from a normal form A, but this normal 
form is diagonal only in the case of the spectral type k = when all eigenvalues are real. 
In the other spectral types the real normal form A is non-diagonal. It has the structure 
discussed in the previous sections. 



4 The paradigmatic example: SL(3; ]R)/S0(1, 2; M) ver- 
sus SL(3;M)/SO(3;M) 



In the present section we illustrate the previously presented theory with a simple, yet 
paradigmatic example, that of the cosets 

SL(3;M) 

■^-^ = S0(p,3-p;M) ^'-'^ 
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where p = 1 is the new case with respect to the well known example of p = (see for 
instance [TO]). From the physical point of view the coset Aio,3 provides the description 
of D = 5 pure gravity reduced three dimensions. Correspondingly the coset A^i,3 is just 
related with time-like brane solutions of D = 5 pure gravity. It can also be related to 
several other interesting brane constructions, yet the viewpoint adopted in the present 
paper is mathematically oriented. We just want to illustrate by means of this examples the 
key points of the explained general constructions. In particular we aim at illustrating the 
relation between the choice of spectral type, the assignment of values to the generalized 
hamiltonians in involution and the analytic structure of the constructed integrals. 

Hence discarding the well known case p = for which we refer the reader to [TO] , we 
concentrate on the novel features of the case p = 1. 



4.1 sod, 2; 



decomposition of the SL(3;M) Lie algebra 



The starting point is to consider a basis of SL(3;]R) generators well adapted to the 
Minkowskian coset SL(3, M) /S0(1, 2; M) rather than to the Euchdian coset SL(3, M) /S0(3; M). 
This is easily constructed with the following procedure. First one multiplies all generators 
adapted to the euclidian basis by the Lorentzian metric tensor 



/ 1 \ 

0-10 

-1 y 



(4.2) 



then one redefines the diagonal Cartan generators by subtracting their trace and making 
them in this way traceless. The result is a set of 8 generators of the s[(3; M) Lie algebra 
with the property that the last three close the Lie algebra of so(l,2;M), while the first 
five span a basis of the spin s = 2 representation of the same. Namely we have: 



[J 

[Ja,K 








^abc 'I 1 



(4.3) 



where Eabc is the standard Levi-Civita antisymmetric symbol and Ii{Ja)ij denote the 5x5 
matrices representing the so(l, 2; M) generators in the spin s = 2 case. 
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Explicitly we have 



and 



Ki 



K, 



Ji 



/ ^ 





o\ 





v 









A 0^ 







/o 





V V2 







/ -1 \ 

-10 




/ -1 \ 

-10 



1 



\ 





V 





/ 






V 3 



\ 



1 

V2 



V 



-71 



/ \ 

0-1 
1 



(4.4) 



(4.5) 



As one can easily note, among the coset generators three are non-compact, while two are 
compact. 



With reference to eq.(3.7) this means that in this case we have 



r = 2 . (4.6) 

Explicitly the non-compact coset generators are the two Cartans Ki 2 and the off-diagonal 
generator K4. The two compact coset generators are instead K35. Similarly the three 
subalgebra generators are distributed into the two non-compact ones Ji,3 and the compact 
one J2. The explicit form of the 5x5 matrices representing the generators Jj is the 
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following one: 



R(Ji) = 



R(J3) = 



^0 





-2 





\ 

















-2 


























-1 










-1 


/ 














^ > 

























-1 











-1 








-1 


V3 












; R(J2) = 



/ 





-1 -V3 

y -10 



1 o\ 



1 








and they are such that 

[J„,K,] = R(Ja)ijKj 

4.2 The Lorentzian Lax operator 



(4.7) 
(4.8) 



The next step is the definition of the Lax operator. According to our established conven- 
tions we define it as follows: 



1 ^ 

L(t) = Y^{t)K, + Y2{t)K2 + -=y2^i(t)K 

v2 



Y2(t) J 



(4.9) 



The so(l, 2; M) Lie- algebra- valued connection is then obtained from the Lax operator with 
the standard i?-niatrix rule 



W{t) = L>(t)-L<(t) 

/ -lY,{t) -lY,{t) \ 

= -lY,{t) -lY,{t) 

\-lY,{t) lY,{t) J 



(4.10) 



and duely satisfies the condition 



r)W^{t) + W{t)r) = 



(4.11) 
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The Lax propagation equation is normalized as follows: 



dt 



Lit) + [W{t),Lit)] = 



(4.12) 



Decomposed along the basis of coset generators K, eq.(4.12) yields the following system 
of five differential equations: 



1 /3 



2 V 2 



n(t) 



2V2 



^^ + -Yi{t) = 
2^/2 dt 



Yi{t)Y,{t) 
V2 



+ 



Yi{t)Y,{t) 
V2 



+ 



d 

dt 
d 

dt 
d 
dt 



lY,it)Y,{t) + j^Y.it) = 0. 





= 0, 


Ysit) 


= 0, 


Y,it) 


= 0, 



(4.13) 



This is the first order differential system on the 5-dimensional Poissonian manifold pro- 
vided by the Borel subalgebra of B(st(3)) C s[(3) which, as a consequence of the general 



discussion of section 2.2 is Liouville integrable. On this 5-dimensional manifold there are 
four hamiltonian functions in involution. Of these latter one is just zero, one is a Casimir, 
that labels the 4-dimensional symplectic leaves into which the 5-dimensional manifold fo- 
liates. The remaining two functions are the 2-hamiltonians which guarantee the Liouville 
integrability of each 4-dimensional symplectic leaf. In the next section we consider the 
explicit form of such hamiltonians. 



4.3 The hamiltonian functions in involution 



As recalled above the differential system (4.13) is Liouville integrable since it admits 



four hamiltonian functions in involution, that can be explicitly constructed according to 



formula (2.49). The first of these four hamiltonians is identically zero since it corresponds 



to the trace of the Lax operator, namely to the extra generator of gl(3, M) which is deleted 
in order to step down to s[(3, M). The last, which, instead of being polynomial, is rational, 
corresponds to the advertised Casimir labeling the leaves of the foliation. 



Explicitly from (2.49) we obtain the following result: 



f)2 = f)02 



YM 



\Y2{tf - -YM 



\Y.{tf - \Y,{tf 



(4.14) 
(4.15) 



1)3 = 1) 



Y^{tf Y,{tYY,{t) , Y,{tfY^{t) , Y,{tfY^{t) Y,{tfY^{t) 



03 



3v^ 



+ 



2VQ 



4v^ 



4^6 
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^ 4^3(t)n(t)n(t) , (4.16) 

h - h - ^iW + ^^W ^3(t)y4(t) 

and by explicit calculation we can verify that the functions f)^ {A = 1, . . . , 4) are constant 



along the Toda flow, namely, upon use of eq.s(4.13) it is identically true that 



dti)A = 0. (4.18) 



If we calculate the secular equation for the Lax operator (4.9) we get 

= Det {L{t) - A 1) 

= -X^ + i)2X + i)3 . (4.19) 
Hence parameterizing the Toda flows by means of the values of the hamiltonians, the 



eigenvalues of the Lax operator are given by the three roots of the cubic equation (4.19). 
Using Cardano's formula these three roots are given by 

2/3 

2V3f}2 + v^2 ( VA - 9fi3 ' 

Ai 



2^[)2 + v^(yA-9f33 



62/3^VA-9()3 
2 (3i + V3)i)2 + v^v^ (1 - iv^) I^VA - 9f)3 



As 

222/335/6^3^^ _ 9(^3 

2/3 

) 

A3 = , ' — (4.20) 

222/335/6 _ 91^3 



2 (-3i + 73) f}2 + (1 + i\/3) (VA - 9f)3 



where 



A = -12[)l + 81i)l (4.21) 
is the discriminant of the cubic equation. 



4.4 Normal form of the Lax operator 



The discussion of the normal forms for the coset 3^(^/2^^^) splits in two cases since, from 
the deflnition of the number p we have 



He 



Indeed the embedding (3.3) is trivial in this case 

U = s[(3,M) ^ s[(3,M) , 
U D =so(l,2;M) ^ so(l,3 - 1;M) C s[(3,M) . (4.23) 
The two possible spectral type are characterized by = or /c = 1. 
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4.4.1 Spectral type k = 



In this case we have from eq. (3.32) T 
operator is 

Ao = Ao = 



V 








To 


A2 





I3 and the normal form of the Lax 








(4.24) 



-Ai - A2 y 



where Ai,2 are two real eigenvalues and — Ai — A2 is the third, also real. From the point 



of view of the Lie algebra, the normal form (4.24) is just a linear combination of the two 
Cartan generators Ki 2. Indeed we can write 



A 



1 
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(Ai - A2) Ki 



(Ai + A2) K2 



(4.25) 



A generic element of the S0(1,2;]R) group can be parameterized as a product of three 
elements of the three one-parameter subgroups, namely we can set 

S0(1,2;M) 3 0{v,e,w) = exp[vJi] ■ exp [0 J2] ■ exp[wJ3] = 

cosh(»;) cosh(io) — sin(6) sinh(i;) sinh(ui) — cos(6') sinh(t)) cosh(io) sin(S) sinh{?;) — cosh(v) sinh(to) 
cosh(i') sin(9) sinh(ui) — cosh(«)) sinh(i)) cos(S) cosh(t)) sinh(t)) sinh(«)) — cosh(v) cosh(w) sin(6) 
— cos(6') sinh(w) sin(6) cos(9) cosh(to) 

(4.26) 

where the parameters v, w and 6 can be regarded as the three Euler angles (two hyperbolic 
and one elliptic) that parameterize S0(1, 2; M). Hence in the spectral type k = the initial 
value of the Lax operator at time t = can be written as 



To(0) 



a ■ Ao ■ O,' 



(4.27) 



where 



Co = 0{v,e,w) . (4.28) 

In this way the matrix Lq{0) depends on the five real parameters {Xi, X2,v,9,w} which 
parameterize the initial conditions Yi{0) for the five real fields Yi{t). Indeed the val- 
ues 1^(0) as functions of \2,v,6,w} can be extracted by projecting Lo(0) along the 
orthonormal basis of coset generators Kj. 

4.4.2 Spectral type k = 1 



In this case we have from eq. (3.32) 



r = Ti 



y 1 



(4.29) 
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and the eigenvalues of the Lax operator are given by a pair of complex conjugate eigenval- 
ues \i = X + iy, \2 = X — iy, while the third one is the negative of their sum, A3 = —2 x. 
Hence the normal form of the Lax operator is as follows: 



Ai 



\ 



Ai 



v 



X + \y 
X — iy 




X —y 
y X 
-2x 

\ 





-2x 



(4.30) 



As element of the Lie algebra, rather then being a linear combination of the two Cartan 
generators, in this case, the normal form is a linear combination of one Cartan and one 
of the compact coset generators. Indeed we have 



Ai 



-VQxK2 - V2yK 



(4.31) 



Hence in the spectral type k 
written as 



1 the initial value of the Lax operator at time t = is 



Li(0) 



(4.32) 



the rotation matrix Oq being defined in eq.(4.28). As a result the initial values Yi(0) of 



the 5 fields are now parameterized by the five parameters {x,y,v,6,w}. 

In both spectral types we can calculate the values of the constant hamiltonians i)A as 
functions of the five real parameter set, either {Xi, X2,v,9,w} or {x,y,v,9,w}. This is 
what we do in the next section. 



4.5 Characterization of orbits through the hamiltonians 

It is now instructive to characterize the orbits and hence the normal forms of the Lax op- 
erators through the values of the conserved hamiltonians responsible for the integrability 
of the system. 

As we have seen the two hamiltonians entering the secular equation and hence the 
determination of the eigenvalues are f}2,3- Furthermore, what distinguishes the two spec- 



tral types is the sign of the discriminant A defined in equation (4.21). When A < 
we have three real eigenvalues, while when A > we have a pair of complex conjugate 
eigenvalues and a third real one. The regions in the ()2, i)3 plane corresponding to the two 
spectral types are visualized in Figjl] It is now instructive to evaluate the explicit form of 
the hamiltonians and of the discriminant in the two spectral types. 
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Figure 1: Phase-diagram in the [)2, [)3 plane. In the locus A > we have complex 
conjugate eigenvalues and spectral type k = 1. In the locus A < there are three real 
eigenvalues and spectral type k = 0. The line A = which separates the two regions is a 
singular locus where two eigenvalues coincide and we have an enhancement of symmetry. 
On this locus the normal form admits a one-parameter stability subgroup. The cuspidal 
point f)2 = f)3 = corresponds to nilpotent Lax operators as we discuss in the paper [36] 
submitted to the hep-th arXiv after the first appearance of the present paper. 



4.5.1 The hamiltonians in spectral type k = 



/^From the definition given in eq.s(4.14 4.17), by using the initial Lax operator (4.27) to 
calculate the fields Yi{0) and hence the hamiltonians, we find 

, 



1)4 



A 



+ A2A1 + A2 , 



(4.33) 
(4.34) 

-AiA2(Ai + A2) , (4.35) 
— [2 (sin^sinhf A^ + 2coshtL'(coshi(7sin^sinhf — coshf sinht(;)A2Ai 
-|- sinh w(sin ^ sinhf sinh w — cosh f cosh w)A2)] 

X [2 (sin 9 sinh v cosh^ w — 2 cosh v sinh w cosh w + sin 9 sinh v sinh^ Ai 
-|-((cosh2w + 3) sin^sinhf — 2coshf coshwsinhw)A2]~^ , (4.36) 



-3(Ai -A2)'(2Ai + A2)'(Ai + 2A2)^ 



(4.37) 
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As it is evident from the above exphcit expressions, in the spectral type A; = the 
discriminant is strictly negative and it reaches the value zero only in the case of degenerate 
eigenvalues, namely when any two of the three eigenvalues are equal. 



4.5.2 The hamiltonians in the spectral type k = 1 



From the definition given in eq.s(4.14 4.17), by using the initial Lax operator (4.32) to 



calculate the fields Yi(0) and hence the hamiltonians, we find 



[)4 



0, 

-2x {x^ + y'^) , 

[cosh V (cos 6 sinh 2wy'^ + 2x cosh w(2?/ sin ^ + 3a; cos 9 sinh w)) — 
sinh V ( (3x^ + y^) cos 9 sin 9 cosh^ w + (3x^ + y^) cos 9 sin 9 sinh^ w 
+ (3x^ + 1/^) cos 9 sin 9 — Axy cos 29 sinh w)] x 

2 ( -X cos 9 sin 6' sinh v cosh^ w + ?/ cosh sin 9 cosh w 

3 

-3x cos ^ cosh f sinh w cosh w H — x cos ^ sin ^ sinh v sinh^ w 

2 



(4.38) 
(4.39) 
(4.40) 



+ -X cos ^ sin ^ sinh f + ycos{29) sinh f sinh 

2\2 



H -1 



A = 12y^ {9x^ + y^) 



(4.41) 
(4.42) 



Once again also in this branch of the solution space the sign of the discriminant is definite. 
A is positive definite and it vanishes only for y = 0. This condition however corresponds 
to real degenerate eigenvalues and matches the same condition obtained from the other 
branch with spectral type k = 0. 

5 Examples of explicit solutions 

In this section we illustrate the integration algorithm by considering some example of 
solutions corresponding to the three spectral types: k = 0, k = 1 and degenerate. 

5.1 An example of solution of spectral type k = 

A very simple solution of this spectral type can be obtained fixing the following initial 
data: 



{Ai,A2,A3} — {| 



3 

2 ' 2 ' 



-2} 



(5.1) 
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On 



exp[2 Ji 



/ 1+^ 

2e2 2e 
2e2 2e2 





{v = 2,9 = Q,w = 0). 



(5.2) 



With these data the initial form of the Lax operator is the following: 

\ / ]^ 1 _-!+€■ 







i^o(O) = Oo 



V 



^ 

" 2 









4e4 



e 
4 



4e4 







4e4 

i(4 + ^ + e^) 
-2 



\ 



(5.3) 



As we see Lo(0) has a block diagonal structure 2 + 1. Such a block structure is preserved 
throughout the all flow from t = — oo to t = +oo as we can deduce from the explicit result 
of the integration 



Yiit) 
Y,it) 



1 _ 26^ + e8 + e2* + 2e2*+4 + e^^+s 
" (-1 + 2e4 - es + e^' + 2e2*+4 + e2t+8) 

2e*(-l + 6^) 



-1 + 2e4 - eS + e2* + 2e2«+4 + e2<+8 



0, 
. 



(5.4) 

(5.5) 
(5.6) 

(5.7) 
(5.8) 



The vanishing of both l4(t) and Y^lt) is what guarantees the block diagonal structure of 
the Lax operator. The same fact however implies that the 4-th hamiltonian, the rational 
one is indeterminate in this case, being the ratio of two zeros. The other two (polynomial) 
hamiltonians have instead the following explicit values: 



13 



3 
'2 



(5.9) 



The plot of the two non-trivial functions 11,3 (t) is exhibited in Fig|2| As we see there is 
a singularity in both fields at a finite time t = to — —0.03663. This singularity separates 
the range of the variable t in two parts. We can consider the solution only on one side 
of the singularity. Let us consider for instance the Cartan field Yi{t). As we know this 
function is actually the derivative of the corresponding Cartan field hi (t) pi5| and in order 
to reconstruct the physical interpretation of our solution we are supposed to perform a 
second integration 



hiit) 



Yi{t) dt 



t - log 1-1 + 2e^ - + e^* + 2e2*+4 + e 
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2t+8\ 



(5.10) 
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Plot of Yi 




Figure 2: Plot of the 11,3 (t) fields in the solution of spectral type k = characterized by 
the following initial data: Ai = |, A2 = |, v = 2, 6* = w = 0. In both fields there is a 
singularity aX t = tg — —0.03663. The other three fields are constant or even zero. 



The singularity at a finite time is a qualitative difference between the type of solutions 
encountered in the pseudo-Riemannian case and those encountered in the case of cosmic 
billiards (Riemannian coset manifolds). A similarity, instead, which exists between the 
spectral type = of pseudo-Riemannian system and the billiard case is the asymptotic 
behavior at t = ±00. For this spectral type (but not for the other), just as in the billiard 
case, the Lax operator tends asymptotically to a diagonal form which differs from A only 
by a permutation of the eigenvalues. We can verify this statement in the present example. 
We find 



Lo(-cx)) 



3 
2 












1 
-2 



7^ A 



35 



Lo(oo) 



/| 

I 



v 







A . 



(5.11) 



5.2 Another solution of the spectral type k = with finite hamil- 
tonians 

To appreciate the differences we apply the integration algorithm to the case where the 
choice of the eigenvalues and of the spectral type remains the same as in the previous 
example but we modify the initial rotation element Oq by switching on also a compact 
rotation angle — j. So we set 



Co = exp [2 Ji] • exp 



TT 



2e2 



2e2 



v 



2V2e2 2V2e2 

2V^e2 2V^e2 

1 1 

a/2 V2 



\ 



(5.12) 



and for the initial Lax operator we get 

/ 3+26^^+36^ 

Lo(0) = OoA^o"' = 



16e4 
3(-l+e'^) 



v 



16e4 



3(-l+e8) 

16e4 

3-2e"'+3e'* 
We* 
7(l+e*) 



_7(-l+e4) \ 

8e2 
7(l+e4) 

8e2 



8e2 



8e2 



(5.13) 



/ 



Calculating the hamiltonians from the above form of the initial Lax operator we obtain 



{f)l,f)2,^3,f)4} 



13 3 1 
'T'~2'~2 



(5.14) 



As wc sec P)2 and [)3, which depend only on the eigenvalues arc the same as before. On 
the other hand, i)4 is no longer undefined as in the previous case and obtains the finite 
rational value — |. This is so because the new initial value of the Lax operator as no 
degenerate minors has in the previous case. 

The new explicit solution is given by the following functions: 



Y3{t) 



(-2 + be"^*) (1 - 2e* + e*^ + 2e2' + e^* + 462*+"* + 2e 



„2t+8 . 



2V2 (1 + e^t) (1 - 2e4 + - 2e2' + , 



3 (_4 + 3e7t) 
2(l + e7«) 



. 4e2t+4 _ 2e2*+8 _ 2e'^*+* + e'^*+®) 



V^e*(-l + e8) (-2 + 5e''*) 



Vl + e^' (1 - 2e4 + - 2e2t + e^"* - 4e2t+4 - 2e2*+8 - 2e''t+4 + e^t+S) ' 

7V2e-«-2 (1 + e4) 

^_e-iit-4 (1 + e''*)^ (1 - 2e4 + - 2e2* + e^* - 4e2*+4 - 2e2*+8 - 2e''*+* + e'^*+») ' 

7(-l + e4) 

'e2^_e-7t-4 (1 + e7t) (1 _ 26* + - 2e2t + e''* - 4e2*+4 - 2e2t+8 - 26^+4 + e'^'+S) ' 



(5.15) 

(5.16) 
(5.17) 
(5.18) 

(5.19) 
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Figure 3: Plot of the Cartan fields Yi^2{t) in the solution of spectral type k = charac- 
terized by the following initial data: Xi = \, X2 = ^, v = 2, 9 = ti/A, w = 0. In Yi(t) 
there are two singularities at finite times t = ts^ — —0.338507 and — 0.0417536. 

The behavior of the solution is qualitatively similar to that discussed in the previous case. 
This is evident from the plots of the Cartan fields exhibited in Fig|3| As one realizes there 
are just two singularities at finite time t = t^i — —0.338507 and tg^ ~ 0.0417536, that 
affect one of the two Cartans but not the other. The same singularities appear in the 
field 1^3 (t). The other two fields Y^j^lt) are real only in the interval between the two 
singularities as it evident from the plot of, for instance, ¥^{1), exhibited in Fig|4j This 
means that the overall solution is properly defined only in the interval between the two 
singularities. Notwithstanding this fact if we calculate the asymptotic limit of the Lax 
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0.05 




Figure 4: Plot of the non-compact nilpotent fields Y4{t) in the solution of spectral type 
k = characterized by the following initial data: Ai = ^, A2 = ^,v = 2,9 = n/4:, w = 0. 
The function l4(t) is real only in the interval comprised between the two singularities 
t = ts,^ -0.338507 and t,^ ~ 0.0417536. 



operator at ±00 we obtain finite diagonal real forms. Indeed we find 



/ 3 

2 



-00 







Lo(oo) 





^ 
^00-2 

} -2 

I 













^ A 



7^ A 



(5.20) 



In the present example the difference between the order of eigenvalues at +00 and at 
—00 is provided by the permutation of highest order, just as it happens in the billiard 
Riemannian case for flows not touching singular surfaces. The previous case did not have 
this property because it developed on a singular surface and the indeterminacy of the 
fourth hamiltonian was a sign of that. 

Although similar to the billiard case the asymptotic limits loose their meaningfulness 
in the pseudo-Riemannian case since they are separated from the physical flow region by 
regions where the Lax operator becomes complex. The real asymptotic diagonal limits 
are approached through imaginary values. The physical flow region is typically bounded 
by singularities. 
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5.3 An example of spectral type k = 1 

As an example of the other spectral type we choose the solution generated by the following 
very simple initial data: 



k 

V 

e 

w 



1 X 

, 

TT 

3 ' 
. 



y 



The corresponding rotation matrix is 



On 



/ 1 




V 







1 _ 

2 

1 

2 2 



V3 
2 



and the resulting initial Lax operator is 

(l 

Lm = 



_1 \ 

2 2 

1 _5 3^/3 

2 4 4 

. V3 3Vl l I 

\ 2 4 4 / 



The resulting vector of hamiltonians is 

{f)l,[)2,[)3,f34} = {0,2,-4,2} 

and the explicit form of the solution is given by 

sec(2t) (Qe^* cos(4t) + 4sin(2t) + 3e6*(sin(4t) + 3)) 



Yi{t) 
Y2{t) 



V2(6e6*cos(2t) + 2) 
(36*5* cos(2t) - 36*5* sin(2t) - 2) 

3e6* cos(2t) + 1 ' 
2^26-3* 

cos 2 (2t) v/2e-6* cos(2t) + 3 cos(4t) + 3 ' 
2V6e"2*(3cos(2t) -sin(2t)) 

^3e2* cos(2t) + 6-4^26-6* cos(2t) + 3 cos(4t) + 3 
2^36* 



cos 2 (2t)v/3e2* cos(2t) + e- 



■4t 



(5.21) 
(5.22) 

(5.23) 

(5.24) 



(5.25) 



(5.26) 



(5.27) 

(5.28) 

(5.29) 
(5.30) 
(5.31) 
(5.32) 



The structure of this solution can be considered analyzing the plots of the various fields. 
The Cartan fields exhibit a quasi periodic behavior (with singularities) displayed in figjs] 
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The two fields 1^3,5 (t) have instead a periodic real behavior for t > tm and for t < tp 
respectively, where 

tp ~ 0.785398 ; tm ^ -0.785398 = -tp (5.33) 

are finite times. Respectively below and above these singularity barriers the fields l3,5(t) 
become imaginary. This is evident from the plots displayed in Fig. [6] This behavior re- 
stricts the physical range of the solution to the interval [—tp , tp]. This is further confirmed 
by the plot of the field Y4{t) which is real only in the same interval. This is seen in FigjT] 

In this way we come to the conclusion that, notwithstanding the appearance of periodic 
functions and the periodic behavior of some of the fields of the system, also in the case of 
the spectral type k = 1, the generic form of the real solution appears to be the evolution 
on a finite range of the time-line, bounded at the extrema by singularities. A similar 
generic behavior is suitable for the description of an evolution from spatial infinity to a 
horizon as it happens in black hole physics. 

A detailed study of the solution space, a classification of the asymptotic limits and 
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Figure 6: Plot of the non Cartan fields 5(t) in the solution of spectral type k = 1, with 
parameters x = l,y = l,v = 0,9 = = 0. 

the analysis of critical surfaces in the moduli space is postponed to future publications, 
where the physical interpretation of the Lax equation solutions in connection with p-brane 
physics will be addressed. 

6 Conclusions 

In this paper we have presented a new view-point on the integrability of supergravity 
cosmic billiards and black holes that is based on the Poissonian structure of the underlying 
solvable Lie algebra S. 

The main results of our paper are two: 

• The explicit construction of the integration algorithm extended also to the case of 
Lorentzian cosets U/H*. 

• The explicit construction of the hamiltonian functions in involution responsible 
for Liouville integrability. 
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Figure 7: Plot of the non Cartan field Y^lt) in the solution of spectral type k = 1, with 
parameters x = l,y = l,v = 0,9 = = 0. Outside of the plotted range the field is 
imaginary. 

We believe that a systematic use of our techniques for the construction of black-hole and 
billiard solutions will provide new results and new insight. In particular the relation of 
the Hamiltonians and the Casimirs with the physical invariants of the solution, like the 
entropy or the total mass will prove very helpful and inspiring. We leave this to future 
coming publications. 

A point which we have not yet addressed but which is of the highest relevance concerns 
the issue of global topology of the solution space. The solvable parametrization covers only 
open branches of this space and the question of how to glue together different branches 
is very important. 

Also this issue is left over for future publications. 
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Note added in the revised version As stated in the previous aknowledgements, 
the authors of [35j had indeed independently adapted Kodama integration algorithm to 
the case of G/H* Lax equations and had derived some explicit solutions for the cases of 
SL(2,M)/S0(1,1) and SL(3, M)/S0(1, 2). The construction of the integration algorithm 
for the case of nilpotent initial Lax operators, which is relevant for extremal Black Holes, 
was performed in full generality in our paper [36j. In a revised version of their paper, which 
appeared the same day as our , the authors of presented some particular solutions 
corresponding to specific nilpotent Lax operators pertaining to SL(3, ]R)/S0(1, 2). 
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